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In the present paper we generalised a phenomenological model developed irA for the description 
of magnetostructural phase transitions and related peculiarities of elastic properties in solid oxygen 
under high pressure and/or low (below 40 K) temperature. We show that variation of all the lattice 
parameters in the vicinity of a/3-phase transition is due to both the shift of basal closed packed 
planes and appearance of the long-range magnetic order. Competition between these two factors 
from one side and lattice compression below T a p from another produces non monotonic temperature 
dependence of lattice parameter b (along monoclinic axis). Steep decrease of the sound velocities 
in the vicinity of T a p can be explained by the softening of the lattice with respect to shift of the 
close-packed planes (described by the constant K2) prior to phase transition point. We anticipate 
an analogous softening of sound velocities in the vicinity of a<5-phase transition and non monotonic 
pressure dependence of sound velocities in a-phase. 
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I. INTRODUCTION 

Solid oxygen that belongs to a family of cryocrystals is being studied for more than 100 years and still attracts 
attention of researches. This simple molecular crystal has very complicated temperature-pressure (T — P)phase 
diagram consisting of magnetic and nonmagnetic, metallic and dielectric, studied and non-discovered phases (see, e.g., 
review^ and references therein). Phase transitions and drastic change of the magnetic, electronic, elastic properties of 
solid O2 can be triggered by variation of temperature, pressure, external magnetic field, etc. The thorough analysis 
of corresponding regularities is rather non-trivial problem because the O2 crystal lattice is mainly hold by the weak 
Van-der-Waals and so is very soft in comparison with ordinary solids. Besides, in contrast to other crystals, exchange 
magnetic interactions at low temperature prove to be of the same order as lattice energy. As a result, application 
of external fields (temperature, stress, magnetic) gives rise to a pronounce variation of all the lattice and magnetic 
parameters. In this situation a precise microscopic calculations should be combined with general thermodynamic 
treatment. 

In the I-st part of this paper— we made an attempt to develop a phenomenologic (Landau-type) model aimed at 
the description of the lattice and magnetic properties of solid oxygen in a sequence of temperature/pressure induced 
[3 — > a — > 5 transitions. 

The proposed model was based on the following assumptions: 

1. Magnetoelastic coupling is so strong that abrupt change of the magnetic structure leads to the noticeable 
variation of the crystal lattice. 

2. The primary order parameter in the series of (3 — > a — > S phase transitions is a homogeneous shift of closed- 
packed planes. The order parameter is defined with respect to the virtual hexagonal, Dg^ pra-phase. 

3. The only macroscopic parameter that controls these transitions is the specific volume. Temperature dependence 
of the specific volume cannot be calculated within the model and should be taken from the experiment. 

As it was shown this model gives good agreement with experimental data. However, some questions have been left 
beyond its scope. In particular, we have considered only two factors that define crystal structure of different phases - 
shift of the closed packed planes and specific volume of the crystal. More thorough analysis should account for rather 
strong rhombic deformations within the plane and separate contribution of interplane distance and isotropic in-plane 
strain into the change of specific volume. 

In the present paper we try to refine the model by taking into account all the parameters that define the crystal 
structure of /?-, a- and 5- phases and explain observed peculiarities of temperature/pressure dependence of lattice 
constants and sound velocities. 



2 



We would also like to mention with our great pleasure that the preliminary results of this paper were presented at 
International Conference CC-2006 (in Kharkov) devoted to the prominent Ukrainian scientist and wonderful woman 
Antonina Fedorovna Prihot'ko who is well known for her brilliant experiments in physics of cryocrystals. She is also 
famous for fundamental results in optical investigation of a and /3-phases in solid oxygen. With the present paper we 
try to pay our tribute of respect to her memory. 

II. MODEL 

It was already mentioned that the crystal structure of /3-, a- and (5-phases can be considered as a homogeneously 
deformed hexagonal lattice. In general, such a deformation can be consistently described with the use of four inde- 
pendent variables: lattice parameters a, b, c, and angle (3 of monoclinic cell. Corresponding combinations that form 
representation of the Dg^ space symmetry group (parent pra-phase) are: 

i) relative shift of closed-packed planes £ = (c/a) cos/3; 

ii) isotropic strain of basal plane As/so = (ab — eto&o)/ ( a o&o); 

hi) relative extension/contraction u zz = (h— ho) /ho in the direction perpendicular to closed-packed planes (Z-axis) 
, where h = c • sin /3; 

iv) shear strain in the basal plane u = (a — v / 3&)/(v / 3a 6o) 1 ^ 2 - 

The quantities with subscript "0" are attributed to a certain reference state (at zero T or P) . In assumption of small 
variation of interplane distance ( formally expressed by inequality u zz <C 1) relative change of the specific volume 
Av/v can be readily expressed in a form of a simple sum Av/v = As/so + u zz . 

Equilibrium values of £, u zz , u, and As/sq at a given temperature T and hydrostatic pressure P are calculated from 
standard conditions for minimum of Gibbs' potential $ that is supposed to be invariant with respect to operations of 
the symmetry group ^>\ h - 

Using results of symmetry analysis given in^ (see Table 2 therein) general expression for <!> can be represented as a 
sum of structural, $ s t r , magnetic, $ m ag, elastic, $ e i a st contributions and interaction term $; n t: 

$ = $ str + $ mag + $clast + $int- (1) 

The structure of the first term in (JTJ 

$ str = -if 2 (s, h) [cos27r£i + cos27r£ 2 + cos27r(£i - f 2 )] 
+iif 4 [cos4<i + cos47r£ 2 + cos4tt(£i - &)] + P f ^ + u zz ^j (2) 

reflects the translational invariance of the crystal with respect to relative shift u of neighboring close-packed planes of 
hexagonal lattice, namely: 27r£i i 2 = bi^u, where bi^ are vectors of the reciprocal pra-phase lattice. Phenomenological 
constants K2 and K4 can be considered as coefficients in Fourrier series of lattice potential. 

The coefficient K2 is supposed to be a linear function of the isotropic in-plane strain and interplane distance with 
corresponding phenomenological constants A s , oc (where A„ = 10 GPa was estimated ini): 

As 

K 2 (v) = K Q - A s XhU zz . 

s 

Both constants A s , A/,. > are supposed to have a positive sign starting from general considerations (confirmed with 
further calculations, see below). Namely, the effective constant K% describes the strength of intermolecular forces 
that keep relative arrangement of basal planes. Increase of the average intermolecular distance should give rise to 
weakening of intermolecular bonds and hence to a decrease of K.%. 
Magnetic contribution in (JTJ 

<f mag = - j(k 13 ) jyp ? - j(k 12) ]T(i^) 2 - j( kl4 ) ]r>f ) 2 , (3) 

i=i j'=i i=i 

accounts for the exchange interaction only, J(kj) > are Fourier components of exchange integrals labeled according 
to stars kj (j = 12, 13, 14) of irreducible representations of Dg^ space group, antiferromagnetic (AFM) vectors 
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\( s ) unambiguously describe the magnetic ordering in (3-, a- and ^-phases 17 . It should be also stressed that 
the magnetic ordering in /3-phase has short-range nature (so called correlation ordering with three sublattice 120° 
Loktev structure^" 7 -) . Collinear long-range magnetic ordering characterized by one of 1W vectors, is established 
in a-phase. An effective Neel temperature TV is close to the temperature T a p of a(3- transition (T a p—23.5 K and 
TV = 40 K at ambient pressure) 2 -, so, saturation magnetization Mo and, correspondingly, |F a )| = Mo noticeably 
depend on temperature in the vicinity of T a p. 
The elastic contribution in ([T]) 



$elas(w) 



1 



'mi 



) 2 + 4uL] + 7^33^ 



(4) 



does not include strain tensor components u xz , u yz , that can be considered as small excitations over finite deformations 
given by the order parameter £\ t 2- An expression ((?]) is valid in the limit of infinitely small strains u xx — u yy , u zz . In 
such an approximation operations of plane shift and deformation commutes, hence, all the parameters of the model 
can be refer to the same "zero" state. Otherwise, one should specify the succession of lattice transformations, use 
different reference frames at each stage and work within nonlinear elasticity approach. 

Interaction energy includes terms that describe different cross-over effects. Here we are concentrated on several 
of them. First, variation of interatomic distances give rise to the change of intra- (characterized with interaction 
constants Aj ntra , Aj_), and inter- (constants A; ntor , Am) plane exchange integrals^. This effect has magnetoelastic 
origin and can be formally described by two expressions. The first one is responsible for isotropic 
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and another one for anisotropic effects 



$ (an) , = -A 

mag— cl 



(a) 



(l( Q) ) 2 cos2< 1 + (1 2 q) ) 2 cos2^ 2 + (4 Q) ) 2 cos27r(a -6) 



A 



A 



A 



(if) 2 cos 2^1 + (l 2 °Tcos2< 2 + (l^ J ) 2 cos27r(a -6 



(' 5 )^2, 



L tt) {{u xx u yv ){\^f + 2u xy [(if) 2 (if) 2 ] } 
f {(u xx - ^)(lf ) 2 + 2u xy [(if) 2 - (if) 2 ] } • 



(6) 



Another coupling effect originates from strong anharmonocity which is a peculiar feature of molecular crystals. Cor- 
responding contribution into thermodynamic potential has a form 



$anhar = ~A e {(u xx - U yy ) COS 2^1 

+ 2 u xy [cos 27T^ 2 - cos27r(£i - £2)]} — A. 



Af 



(7) 



where the coupling constants A e , X± are of the same nature as A s , but from the general point of view should be 
much smaller in value (see Table HJ. Really, \ s , Ah describe variation of average interatomic distances that result 
from reconstruction of crystal lattice in the course of phase transition. These constants are proportional to Grunaisen 
parameter which is very high. At the same time A e and Aj_ are responsible for a "differential" effect, namely, variation 
of anisotropy of the crystal lattice in the course of phase transition. 

Standard minimization procedure enables to obtain the value of lattice parameters and stability conditions of (3— , 
a— and <5— phases. 



III. TEMPERATURE/PRESSURE DEPENDENCE OF LATTICE PARAMETERS 

Equations for order parameter and stability regions were derived and discussed in details ini. The refinement of the 
model (splitting of in- and out-of-plane contribution into specific volume) gives rise to no qualitative changes. In this 
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Table I: Phenomenological coefficients 



Coefficient 


Value 


Comments 


Xs 


0.31 • 10 cm /dyne 


at amb. P, 




3.2 ■ 10 cm /dyne 


at 1 =19 K, ^ = 1 -j- 10 GPa 


I3p(T) 


1.6 • 10 1/K 


Tj23.5 K 




i t 1 n — 3 i /TV 

1.7 • 10 1/K 


26.0 Kji 


K cS 


3 GPa 


calc. at amb. P 


K4, 


5 GPa 


calc. at amb. P 


As 


lo bra 


at amb. _P, 




no fPn 


p+ T— 1 Q PC P — 1 -in PPfl 

dl _f — XV, / — _L — _LU \jL (4 




10 GPa 


in ncsi i tyi Tit i nn i~ n /^i o — 1 f) rlPa 

111 dM Ulll U L1U11 UlldL L j j _1 L/ v. I 1 CI. 


Ax = C13 


1 GPa 


in assumption that 033 = 10 GPa 


A e 


0.06 GPa 


in assumption that c'=l GPa 


A^M 2 


-0.04 GPa 




A^Mq 2 


-0.02 GPa 




Aintor Mq 


oc 0.06 GPa 


in all phases 



section we discuss only some additional effects that can be explained and predicted in the framework of the model. 
As ini we consider a homogeneous (from crystallographic point of view) phase with £ = £1 = £2/2. An effective 
macroscopic order parameter that vanishes in /3-phase is introduced as 77 = 1 + 2 cos 27r£. Magnetic ordering in a— 
and 5- phase is described by a single AFM vector and 1^, correspondingly. 



A. Order parameter and isomorphic in-plane strain 



In the "more symmetrical" (3- and S- phases an order parameter takes the limiting values r/ — (/3-phase) and 
77 = — 1 ((5-phase). Saturation magnetization Mq (or correlation parameter in /3-phase) are constant. Pressure and 
temperature dependence of the isomorphic in-plane strain As/s in these cases cannot be calculated from general 
thermodynamic considerations and is determined experimentally as 



— = f T MT';P)dT'- Xs (T)P+[A[^ 
s Jo 



(8) 



/3-ph (5-ph 



where the in-plane thermal expansion coefficient (3 S (T; P) and isothermal Xs(T) compressibility vary in a wide range 
with pressure and temperature (for example, Abramson et al observed 10% change of compressibility per 1 GPa in 
/3-phase at room temperature^). 

In a-phase the temperature/pressure dependence of rj and As/s can be calculated from the system 



K ~K 4 + Aj Q) M 2 - \ s — - \ h u z 



As 



V + 2^(3r, 2 



A[, a) M 2 =0, 



— = C AC*"; P)dT> Xs P +^1(^-3)- Xs A^ a M 2 (T). 
s Jo * 



(9) 



(The behavior of interplane deformation u zz will be discussed in the next section) . The values of phenomenological 
constants derived from fitting of functions 77 (T, P), As/s(T, P) to experimental data^ifi with due account of Mq(T) 

temperature dependenc e 11 ' 12 are given in Table |TJ It can be seen from the Table [J that K eS = K - K 4 + A < f ) M§ and 

K4 are comparable in order of value with the shear modulus of the material^. The value of Xs which we associate 
with the "seed" compressibility is much smaller compared with experimentally observed values 2 -. This means that the 
main contribution into compressibility arises from shifts of close-packed planes and magnetic interactions. 
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B. Interplane distance 



Distance h between the close-packed planes can be calculated from equations 



A i As A 



C33 s C33 

A± As _ X h 
C33 s 2c 33 
A_L As X h 
C33 s 2c 33 



lnter M 2 for 0- phase, 



A 



(a) 



V 



C33 

^£tcr M 2 
C33 



-Mq(T) for a -phase, 



A 



for <5 — phase, 



(10) 



once the dependencies rj(T,P) and As/s(T, P) are known. Analysis of equations (|T0|) shows that in the /3-phase 
temperature/pressure dependence of u zz (and hence interplane distance) is related with lattice compression within 
the close-packed plane. Thus, constant A^ coincides with the elastic modulus C13. Estimated ratios C13/C33 = 0.1 
at ambient pressure and C13/C33 = 0.24 at T — 19 K, P = 1 -j- 7 GPa and also the positive sign of C13 > are in 
agreement with observations made by Abramson et alA 

Figure [T] shows temperature dependence of interplane distance experimentally measured (squares) and calculated 
from (|10p (fitting parameters calculated in assumption that C33 = 10 GPa are given in Table U). Below T a p an 
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Figure 1: Temperature dependence of interplane distance (squares) calculated according to experimental data^. Theoretical 
curve (solid line) is calculated from Eq. (|lUp . Arrow indicates the point of a/3- transition. 



interplane distance abruptly diminishes and then its temperature derivative changes sign. This fact may be explained 
by the combined influence of the magnetic ordering (characterized by difference A^? er — Afc^) and shift of basal 
planes (corresponding interaction constant Xh)- In /3-phase each O2 molecule is situated over the center of underlying 
triangle. In-plane compression (that arises due to cooling or application of pressure) pushes out the overlaying 
molecules in vertical direction, so interplane distance increases. In a-phase hexagonal planes are not only shifted from 
an ideal "close-packed" (with respect to underlying layer) position but can also "slip" under compression. Appearance 
of an additional degree of freedom allows lattice compression not only in-plane but also in perpendicular direction. 
Formally this means that constant Xh is positive, in agreement with monotonic decrease of interplane distance under 
pressure- induced compression (see Fig. [5]). 

Contribution of magnetic interactions into interplane distance is important not only at a/3-, but also at ad- transition 
as can be seen from Fig. [5J If we assume that a- and 6- phases have the same magnetic structure, then pressure 
dependence of interplane distance according to (|10p should be continuous up to (hypothetical) 2-nd order transition 
point at P2 — 7.5 GPa (dashed line in Fig. [J). Change of magnetic structure that originates from inter- plane exchange 
interactions induces I-st order phase transition at Pi = 6 GPa< P2 followed by abrupt change of the order parameter 
77 and all the related parameters including u zz . So, difference between the observed and hypothetical u zz value above 
Pi (see Fig. [5] and Eq. (fTU|l ) is proportional to A[^ C1 . 



Aj" t ' er and is of magnetic nature. 
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Figure 2: Pressure dependence of interplane distance at T=19 K: (squares) calculated according to experimental datai^. 
Theoretical curve (solid line) is calculated from Eq. f|10p . Arrow indicates the point of the real, Pi, and hypothetical, P2, a8- 
transition. 

C. In-plane deformation: parameters a, b 

It is quite obvious from symmetry considerations that establishment of long-range magnetic order in a-phase is 



followed by in-plane deformation described by u x 



^yy 



u. Analysis of the expressions ©, ® shows that the same 



effect can be produced by the shift of basal planes (term with A e ) as seen from the following equation: 



-j(v -3??) + • 



(11) 



Pressure dependence of u xc 
-0,11 



calculated from Eq. (fTTj) along with experimental data is given in Fig. [3] It 



7=19 K 




P, GPa 



Figure 3: Pressure dependence of in-plane rhombic strain u x 
data^. 



(solid line) calculated from Eq. (|f f [) . Squares: experimental 



is clearly seen that in a-phase lattice anisotropy u xx — u yy monotonically increases in absolute value with pressure. 
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This effect can be explained from very simple considerations. Pressure produces compression in all directions in the 
basal plane but relative shift of planes can induce an effective tension in b direction as seen from Figd] (O2 molecules 
move apart in b direction when overlaying molecule moves toward the edge). As a result deformation in basal plane is 
essentially anisotropic and anisotropy is more pronounced at smaller interatomic distances (i.e., at higher pressure). 



a 



-o 



Figure 4: In-plane displacements (arrows) of O2 molecules induced by relative shift of the overlaying (open circle) plane. 

The same effect takes place during cooling a-phase at ambient pressure. Fig. [5] shows temperature dependence of 
a and b lattice parameters below T a p. Intcrmolecular lattice distance a decreases monotonically because of cooling- 
induced compression. In contrast, temperature dependence of the parameter b is non-monotonic. Increase of b 
during cooling from 23.5 to 18 K means that tension in this direction is stronger than cooling-induced compression. 
In particular, Eq. (jlip demonstrates that two mechanisms may be responsible for this behaviour: repulsion of 
ferromagnetically ordered neighbors^ (described by the term A±M$(T), Mq increases with temperature decrease) 
and already mentioned shift of basal planes. 



3,427 




T, K 



Figure 5: Temperature dependence of lattice parameters a and b (solid lines) calculated from Eq., (JTTJ. Squares: experi- 
mental data£). 



Comparison of the above theoretical dependences with experimental data obtained by different groups using different 
technique makes it possible to estimate the values and range of phenomenological constants (see Table fl} . All the 
coefficients could be grouped in three categories: i) large (oc 10 GPa) constants responsible for anharmonicity effects 
{\ s ,h)', h) intermediate (oc 1 GPa) constants (K c s, K4, C13) that characterize elastic properties of crystal and iii) 
small (below 0.1 GPa) constants that describe magnetic interactions and anisotropy effects (rest of the constants). 
Though small in value the magnetic interactions reveal themselves in the case of competition between different 
structural interactions like pressure induced compression and repulsion that arises due the shift of basal planes. 
Strong anharmonicity is quite natural phenomena for molecular crystals. It is remarkable that at ambient pressure 
both coefficients A s and Xh are of the same order and the main contribution into effective structural constant K2 
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arises from isomorphic in-plane strain^ As/s. 



IV. PECULIARITIES OF SOUND VELOCITIES IN THE VICINITY OF PHASE TRANSITION POINTS 



Experimental study of sound velocity in solid oxygen provides information about interactions that play the leading 
role in phase transitions. High pressure measurements in /3-phase 8 show monotonic increase with pressure of all the 
elastic constants except shear modulus C44, the value of which abruptly decreases in the vicinity of /^-transition point 
(approximately 8 GPa at 295 K). Low-temperature curve o 13 i 14 have a pronounced minima at T a p. 

Critical behaviour of sound velocities can be quite naturally explained in the framework of the developed model 
using the concept of Goldstone mode. Reallythe Il-nd (or weak I-s1j2S) order phase transitions are usually accompanied 
with softening of a certain bond responsible for appearance of the order parameter. So, characteristic frequency of the 
corresponding excitations (symmetry related to the order parameter) should vanish or at least noticeably diminish in 
the vicinity of the phase transition point (Goldstone mode). In the case of (3 — > a — ► 8 transitions an order parameter 
(shift of basal planes) is symmetry related with strain tensor component u zx (see Table 2 ini) and, as a result, with 
the transverse acoustic modes propagating in [0001] and [1000] directions. Thus, these modes should have peculiarity 
in phase transition points. 

Sound velocities of the "soft" transverse, v tl and longitudinal, vi acoustic waves propagating in [1000] direction are 
expressed through the elastic modula in a standard way 



1 

1 

2i> 



en + c 55 - J (cix - c 55 ) 2 + 4c 2 5 



en + c 55 + ^(cn - c 55 ) 2 + 4c 2 5 



(12) 



where p is crystal density. 

Elastic modula that equal to the 2-nd derivatives of thermodynamic potential $ (see Eq.JT])) with respect to strain 
tensor components are calculated at equilibrium values of lattice parameters: 



en 



<9 2 $ 



d(As/s 



Cl5 



csin/3 <9 2 $ 



a d^d(As/s) 



du„ 



csin/A 2 <9 2 $ 
a ) ~dl? 
csin/3 



(13) 



Factor (c/a) sin /3 in (fL3"|) is deduced from geometrical considerations. 

Then, omitting some terms immaterial for further discussion we obtain that in all three phases en = xj 1 , c i5— 
in (3- and <5-phases and 



C55 



4nc sin/3 



\ S ((3 S T - XsP) - ^^eff 



(14) 



in 13- phase, 



C55 



47rc sin /3 



X s (XsP - PsT) + K cS -Ki + 2A[, Q) Af 2 



(15) 



in 5-phase, and 



C55 



CIS 



/ 47rcsin/3 

V a 
47rcsin [3 



X s ( Xs P - PsT) - K cS - 3K iV (l - |) 



X s sin 27r^ 



sin 2 2?r^, 



(16) 



in a-phase. 

It can be easily seen from (fl"4f that in /?-phase shear modulus increases with temperature and decreases with 
pressure, at least in the vicinity of transition point, in consistency with observations 8 . 

More interesting is comparison of experimental and theoretical temperature dependencies 21 c tj / = pv\ l shown in 
FigsjSl [7] Anomalous softening of ct agrees well with theoretical predictions and thus may be explained by high 
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T,K 

Figure 6: Temperature dependence of shear modulus c t calculated from (|12|l - (|16l ) (solid line). Squares: experimental data^. 



5,2-, 




10 15 20 25 30 35 40 



T, K 

Figure 7: Temperature dependence of compression modulus ci calculated from (|12|l - (|16|) (solid line). Squares: experimental 



compliance of O2 crystal lattice with respect to shift of the closed-packed molecular planes. Softening of q is not so 
obvious from general point of view but in the case of molecular crystal may originate from strong anharmonicity. 

An analogous softening of elastic modula and corresponding sound velocities is also expected in the vicinity of 
a<5-transition. Fig. [8] shows hypotetical pressure dependence of the modula product^ aci calculated on the basis 
of equations p^)l - (fTB|) with characteristic values of phenomenologic parameters (see Table [XJ) . It is clear that this 
dependence should be nonmonotonic, with noticeable decrease of modulus while approaching to a/3 and a<5-phase 
boundaries. 
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Figure 8: Pressure dependence of the product of CtCi calculated from (|12[1 - (|16|) (solid line). 



V. CONCLUSIONS 



In summary, we have calculated temperature and pressure dependence of all the crystal lattice parameters of solid 
oxygen in the magnetic a/36 phases. The above phenomenological model ascribes the leading role in structural changes 
to the shift of satisfactory agreement with the experiment. 

In particular, shift of the planes is responsible for nonmonotonic temperature dependence of lattice parameter b in 
a-phase (effective repulsion), change of the effective thermal coefficient in c direction in af3 transition (nonmonotonic 
t dependence of interplane distance), softening of shear modulus and corresponding velocities in the vicinity of T a p. 

We assume that with certain stipulation the developed model may be applied to interpretation of IR spectr a 15 ' 16 
which pressure dependence looks similar to pressure dependence of the order parameter. In particular, we predict 
non-monotonic pressure dependence of elastic modula in a-phase and critical behaviour in aS- transition point. At 
the same time some of the problems are still being outside of theoretical treatment. First of all it concerns magne- 
tocrystalline structure of e-phase, metallization of highly compressed oxygen along with corresponding mechanism of 
superconductivity, optical properties of high pressure phases, etc. These and other problems are the subject of future 
investigations. 

We would like to acknowledge all the participants of CC-2006 Conference for keen interest to our presentations and 
fruitful discussions. 
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